Theory of optical spin orientation in silicon 
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We theoretically investigate the indirect optical injection of carriers and spins in bulk silicon, using 
an empirical pseudopotential description of electron states and an adiabatic bond charge model for 
phonon states. We identify the selection rules, the contribution to the carrier and spin injection in 
each conduction band valley from each phonon branch and each valence band, and the temperature 
dependence of these processes. The transition from the heavy hole band to the lowest conduction 
band dominates the injection due to the large joint density of states. For incident light propagating 
along the [001] direction, the injection rates and the degree of spin polarization of injected electrons 
show strong valley anisotropy. The maximum degree of spin polarization is at the injection edge 
with values 25% at low temperature and 15% at high temperature. 

PACS numbers: 72.25.Fe,78.20.-e 



I. INTRODUCTION 

The optical injection of carriers is a powerful method 
for the study of the properties of semiconductors, and 
the optical injection of spins, i.e., optical orientation, 1 
is an important element in the toolkit of the field of 
spintronics. 2-4 Most research has been focused on direct 
optical transitions. There have been fewer studies of indi- 
rect optical transitions, in which the excited electron and 
hole have different wavevectors, and phonon emission or 
absorption processes are necessary to conserve the total 
wave vector. In this paper, we consider the indirect op- 
tical injection of carriers and spins in bulk silicon. 

Macfarlane et al£ first measured the fine structure of 
the absorption-edge spectrum in intrinsic bulk germa- 
nium and silicon in 1958, and identified the phonon- 
assisted indirect gap absorption branches. Thereafter, 
these studies were extended to doped silicon^— The- 
oretically, Elliottii was the first to study indirect ab- 
sorption of excitons under the effective mass approxima- 
tion, and identified the absorption lineshapes at pho- 
ton energy near and far away from the indirect gap. 
He found that the lineshapes are not sensitive to exci- 
ton effects at high photon energy. The electroabsorp- 
tion in indirect gap semiconducto r 12 ! 13 and the absorp- 
tion spectra of multiexciton-impurity complexes 1 ^ were 
also investigated. Hartmani^ determined the absorption 
spectra using a parabolic band approximation. Dunn 16 
and Cho w 1 ' 18 used a Green function method to inves- 
tigate the physical processes in indirect absorption. All 
these models approximated the matrix elements of the 
electron-phonon interaction by their band edge values. 
Later, pseudopotential models 1 ^— were employed to cal- 
culate the transition matrix elements of the electron- 
phonon interaction around the band edge. However, a 
full band structure calculation of the full spectrum of in- 
direct gap absorption is still absent, even with the neglect 



of the excitonic effects. 

Investigations of indirect optical spin injection are less 
common than those of indirect gap carrier injection. 
They are also less common than those of direct gap spin 
injection^ despite the fact that the first optical orien- 
tation experiment^ performed by Lampel in his study 
of the nuclear polarization of 29 Si in bulk silicon, em- 
ployed indirect absorption. The degree of spin polariza- 
tion in such an injection process is sometimes understood 
as a spin-dependent virtual optical transition combined 
with a spin-independent phonon absorption or emission 
process i 26 ' 27 However, due to spin-orbit coupling the elec- 
tron states are not pure spin eigenstates, and the effect of 
the electron-phonon interaction on the indirect gap injec- 
tion needs to be calculated in detail. While Li and Dcry 28 
have recently studied the degree of circular polarization 
of the luminescence associated with the recombination 
across the indirect band gap following the electrical in- 
jection of spins in silicon, a theoretical investigation of 
indirect optical spin injection is still absent. 

In the present paper, we perform a full band struc- 
ture calculation of the indirect optical injection of carriers 
and spins using an empirical pseudopotential model 2 ^— 
(EPM) for electron states and an adiabatic bond charge 
model'' 2 (ABCM) for phonon states. Compared to k ■ p 
models and ah initio models, which are widely used for 
direct gap carrier and spin injection calculation , 24 ' 33 the 
advantage of the EPM is that the electron-phonon in- 
teraction in the whole Brillouin zone can be calculated 
consistently in combination with the ABCM. This ap- 
proach has been successfully used to describe spin relax- 
ation processes^ 4 - and photoluminescence^ 8 . in bulk sili- 
con. 

We focus on the calculated optical indirect injection 
coefficients of carriers and spins, identifying the contribu- 
tion from each valence band and phonon branch. We take 
the excited electrons and holes as free carriers. Because 
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the electron- hole interaction is nearly spin- independent, 
to good approximation it should affect the spin injec- 
tion and the carrier injection rates in the same way, and 
not affect the degree of spin polarization (DSP) of the in- 
jected electrons, which is the ratio of these two quantities. 
The dependence of the injection coefficients and DSP on 
photon energy, conduction band valley, and temperature 
are established. For the injection of carriers, our numer- 
ical results agree with the experiments at high photon 
energy. In the course of our investigations we also dis- 
cuss in detail a simple but widely used model, in which 
the values of the transition matrix elements are approxi- 
mated by their values at the band edge, and we compare 
its predictions with our calculations. 

We organize the paper as follow: We begin with the 
general formula for indirect carrier and spin injection as- 
sisted by phonon emission and absorption in bulk silicon 
in Sec. HH Then the selection rules that follow from the 
crystal symmetry are discussed in Sec. Mil We intro- 
duce the EPM and ABCM models in Sec. [TY1 Finally we 
present our results and conclusions in Sec.lVland Sec. lVIl 
In an Appendix, we describe an improved adaptive lin- 
ear analytic tetrahedral integration method (LATM) 24 i 35 
that we use to perform the six-fold integration over the 
Brillouin zone (BZ). 



II. MODEL FOR CARRIER AND SPIN 
INJECTION BY INDIRECT ABSORPTION 
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FIG. 1: (color online). Band structure of silicon, calculated 
from the EPM, with symmetry notation at the F point and 
the conduction band edge (Ai). The valence band is noted 
by HH, LH, and SO bands. Two important phonon-assisted 
transitions from valence band edge (red dot) to conduction 
band edge (red filled square) are given by green lines and 
blue lines. The dotted lines refer to optical transitions, while 
the dashed lines refer to the transitions mediated by phonons. 
The red hollow square (circle) stands for one of the possible 
intermediated states in the transition given by blue (green) 
lines. 



Fig. Q] shows the band structure of silicon in the en- 
ergy range [—12 eV, 6 eV] calculated from the EPM; 
the details are given in section IIVI The lowest four 
bands are valence bands, of which the upper three are 
the heavy hole (HH), light hole (LH), and spin split- 
off (SO) bands. The valence band edge (red dot) is 
at the r point with = 0. All other bands are con- 
duction bands. From the figure, the conduction band 
edge (red filled square) is at k® « 0.85rX on the A 
symmetry line, and results in six equivalent valleys that 
can be denoted as X, X, Y, Y, Z, Z, indicating the loca- 
tion of the valley center. The calculated direct band gap 
at the r point is E g = 3.4 eV, while the indirect band 
gap is Eig = 1.17 eV. When the photon energy satisfies 
Ei g < Hlu < Eg, optical injection occurs only across the 
indirect gap. Because the excited electron and hole have 
different wave vectors, the transition must be assisted 
by phonon emission or absorption. The green and blue 
lines show two possible indirect gap transitions between 
the valence and conduction band edges. The red hollow 
square and circle indicate possible intermediate states. 

For an electric field E(i) = E^e -1 ^ + c.c, the carrier 
and spin injection rates in silicon can be generally written 
as 

h(T,n) = C b (T,uj)E«(Eiy , 
S'(T,w) = ( fab (T,u;)E:{Ety . (1) 

Here £ ah (T,w) and (f ab (T,u) are the injection coeffi- 
cients for carriers and spins, respectively, at temperature 
T. The superscript Roman characters indicate Cartesian 
coordinates, and repeated superscripts are to be summed 
over. 

Using Fermi's Golden Rule, we find the injection coef- 
ficient £ ah and £^ ab to be of the form 

A ab = A tvx± , (2) 

T;cv\± 

with 

k c el.k v 

X N {k c -k v )\±A a X cV k v ;\ 7 (3) 

where 

J ^ck a vk v -\ = (c ' k c \A\ck c )T£ kc y kvX {Tffk a vk v \) ■ 

(4) 

Here I is the valley index; c(y) is the conduction (valence) 
band index without including spin; k c (k v ) is the electron 
(hole) wave vector, where k c G / means the summation 
is over the I th valley; e c fc c (e„fc„) give the energy spectra 
of conduction (valence) bands; Ml q \ gives the phonon 
energy at wave vector q and mode A (longitudinal optical 
(LO) and acoustic (LA), and transverse optical (TO) and 
acoustic (TA) branches); and N q \± = JV^ + jij, where 
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N q x is the equilibrium phonon number. The operator A 
in Eq. (J4j) stands for the identity operator / in the carrier 
injection calculation, and the / component of the spin 
operator in the spin injection calculation. The indirect 
optical transition matrix elements are 



rpa _ e v * 

ck c vk v X * 



W — 0J nvk 



Ldcnk c 



(5) 

in which c = {c, cr c }, d = {c,a' c }, v = {v,a v }, and 
n = {n,a n } are full band indexes, with a c ,a' c ,a v be- 
ing the spin indexes; e nk is the electron energy at 
band n and wave vector k, and u} nm k is defined by 
Swnmfc = £nfc — £mk- The velocity matrix elements are 
given by Vnmk = (nk\v\fhk) , with the velocity operator 
v = dH e /dp and the unperturbed electron Hamiltonian 



H e ; M nk 



(nk c \H e x p (k c 



k v )\mk v ) are matrix 



elements of the electron-phonon interaction H ep , which 
is written as H ep = J2 q \ H^{q){a q x + a ] _ qX ), with a qX 
being the phonon annihilation operator for wavevector q 
and mode A. 

We now turn to the symmetry properties of Aj b cv \± ■ 
Though bulk silicon has Oh symmetry, each conduction- 
band valley only has C^ v symmetry .~ Therefore each 
Ai-cv\± only has C$ v symmetry, since the summation 
over k c is limited to k c in the I valley. The summation 
of Eq. ([3]) can be rewritten as 



tab 



^-N v N c j ^2 E TT E s ( s ok c ~ £ vk v 

±^(fc c -P„fc„)A - ^) N (k c -P v k v )\± 
iab 



v A a0 

X ' / V:efc e i;(P„fc„)A > 



where 



lab 
•^I-.ckcVkvX 



cj 



ih c )v(P c ,ik v );X • 



(6) 



(7) 



The prime indicates the summation is only over the ir- 
reducible wedge of the Brillouin zone, P c j are the N c ,i 
symmetry operations in C^ v that keep the I th valley un- 
changed, while P v are the Af v symmetry operations in 
Oh- For each process {J; cwA±}, the symmetry properties 
of the tensor with components Aj b cvX± are the same as 

those of the tensor with components A ab ck vk A . There- 
fore, in the I = Z valley, the £§ b cl , A± form a second rank 
tensor with only two nonzero independent components, 



pxx _ tvy _ Ai) 

t,Z;cv\± ~ SZ;euA± — S CT A± > 



i 



Z;cv\± 



- t(2) 

— ?cuA± 



(8) 



Similarly, the Cz b cv\ f° rrn a third rank pseudotensor with 
only two nonzero independent components, 



>-ZXT) 



Z;cvX± 



/■zyx r (i) 



/ XV z 

c 



Z:cv\± 



f-xzy 

^Z:cv\± 



^yxz 



Z;cv\± 

j-yzx 
'^Z;cvX± 



-iC {2) 



(9) 



where C (2) and C (1) are real numbers, and is 
also a real number because of inversion and time reversal 
symmetry in bulk silicon. The injection coefficients in 
other valleys can be obtained by properly rotating the 
Z valley to the corresponding valley. The total injection 
coefficients 



/tab 
• Al cvX± 



= E^I 



cuA± 



(10) 



have higher symmetry, and the nonzero components sat- 
isfy 



pxx _ pyy 



tzz 
$cv\± 



j-xyz /-yzx /-z%y 

^cuA± ^>cuA± ^cfA± 



2£ (2) 

z ?cuA± 
fXzy 

* (2c£L - 



. u (1) 



fzyx 



4C (2) 



*yxz 



All quantities keep the identified symmetry properties on 
summation of one or several subscripts in {cv\±}. 



III. TRANSITIONS AT THE BAND EDGE 

The values of the matrix elements Tck c vk v ,X at the 
band edge, T ek o^ k o ,\, provide insight into the importance 
of each injection process. For the indirect gap injection 
in silicon, the conduction band edge in the Z valley is 
at fc° = (0,0, /ca), while the valence band edge is at the 
r point, fey = 0. In analyzing transitions near the band 
edge, we can use the following two approximations: i) 
the values of the transition matrix elements [Eq. @] are 
taken at their band edge values, and ii) the values of the 
wave vector, the energy, and the phonon number of the 
phonons involved are all taken at their band edge values. 
Under these approximations, Eq. (JB|) becomes 

At CVT± {T,u)^^J cv {hwThn k « x )N k o T± A% VT , (11) 



with 



j cv (hu) = ^ s ( £ ck c 



~'vk v 



Jab 1 V^V^ lab 

■^I-cvt — j^j- / j / ,- M -ck°v(P„k°. 



Hlu) 



)A 



(12) 
(13) 



A£r P„ 



Here r indicates the phonon branches (TA, TO, LA, LO), 
J cv (hu}) is the joint density of states (JDOS) for indirect 
gap injection, Aj b CVT gives the symmetrized transition 
matrix elements at the band edge, and X^Aer indicates 
summation over all modes in the r th branch. Under the 
parabolic band approximation^ the JDOS is 



J cv (hhj) oc (E ig - frioY 



(14) 



We now turn to A] b CVT . Because the HH and LH 
bands are degenerate at the T point and their wave func- 
tions are not uniquely determined, unambiguous values 
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of A^l 0vk0X and Af k a vk Q X for the HH and LH bands sep- 
arately do not exist at the band edge. But this is not a 
problem for Aj, CVT , due to the summation over all sym- 
metry operations. To show this clearly, we indicate all 
intermediate states in the transition matrix elements im- 
plicitly by writing 



' ck c vk v ; A 



with the operator 



— (ck c \T ckcVkv . x \vk v ) 

U! 





phonon mode A 


TA/TO 


LA 


LO 


X 


y 


z 


2 2 

x -y 


\X) 





(0,0, if') 


(0,T 3 ,0) 


(14,0,0) 


\y) 


(0,0, If) 





(T 3 ,0,0) 


(0,T 4 ,0) 


\z) 


(0,T 2 W ,0) 


(if, 0,0) 





(o,o, r 5 ) 



(15) 

TABLE I: The transition matrix elements (Tf fe o „ fc o A , 
T ck°,vk°,A> T ck°,vk°,\)- The symbols with superscript prime 
are for TO phonons. All quantities can be taken as real num- 
bers. 



£ck c — H e — hw A 



huj - H e + e v k v 

H c x p (k c -k v ) (16) 



(|s,±s)) bands with the following states: 



Substituting Eq. ([15]) into Eq. (@|, it is easy 
to find that the expression for Af! CVT includes a 
summation N V ^ 1 ^ P P v \vk®) {vk®\P v , which equals 
Es'glhhh Wk°)(v'k v \/\ when v is the HH or LH band. 
So A^ mT = ^?f cLHT = (^f cHHr + A% WT ) /2 are unam- 
biguous, and give the same value for both the HH and 
LH bands. This conclusion mirrors a similar one in the 
study of injection across the direct gap^ In the follow- 
ing, we focus on selection rules for A a ^. cvX . Li and Dery 28 
also discussed the selection rules in the context of lumi- 
nescence by considering only the lowest conduction band 
and the highest valence band as the intermediate states; 
they also assume equal amplitudes for the two interfer- 
ence processes shown in Fig. [TJ Here we give a general 
discussion using Eq. (|15[) . without relying on the prop- 
erties of the intermediate states. 

Without spin-orbit coupling, the valence band states 
at the r point transform according to the representation 
-^25 (with basis functions that transform as {yz, zx,xy}, 
which are labelled as {X,y,Z} here)^ We first con- 
sider the electron states that lie in the Z valley, in which 
the conduction band edge state transforms according to 
the representation Ai (basis function {z})£H The phonon 
states involved transform according to the representa- 
tion Ai (basis function {z}) for the LA phonon mode, 
A 2 (basis function {x 2 — y 2 }) for the LO phonon mode 
and A5 (basis functions {x, y}) for the TA/TO phonon 
modes H x p (k®) has the same symmetry properties as 
the X th branch phonon polarization vector. The velocity 
v transforms according to the representation Ai © A 2 . 
All nonzero components of the vector T ck o vk o\ are listed 
in Table U for the different valence bands and phonon 
modes. In Li and Dery's approximations, T 2 = 2T\ and 
T 2 = 2T[ because the two transition processes summed 
are of equal magnitude, and T5 vanishes due to their lim- 
itation of the intermediate states. In contrast, we keep 
all nonzero terms in our treatment, and give their values 
in Sec. |lVl 

With spin-orbit coupling, the valence bands at the T 
point are split into HH (||,±|)), LH (||,±i)) and SO 



1 

'2'' 
1 

'2'' 
3 

'2'' 
3 

'2'' 
3 

'2'' 
3 

'2'' 



±=\ Z )\v + ±=\x + iy)\i) 
^|*-0>>|t)-^)U) 
l\z)\f)-±=\x + iy)\i) 



-L\x-iy)\t) + 

--L|A" + ^>|t) 
±\x-iy)\y . 



\\z)\i) 



(17) 



The SO band at the V point is lower than the degen- 
erate HH and LH bands by 44 meV. At the conduc- 
tion band edge, the states can be approximately writ- 
ten as \z f) and \z J.) due to the very small spin mixing. 
The transition matrix elements between these states with 
spin-orbit coupling can be obtained by linearly combin- 
ing the terms in Table [I] 

For a~ light, E = (x — iy)/y/2, we show in Fig. [2]thc 
possible optical transitions from the valence bands to the 
Z valley of the conduction band. Here the spin quanti- 
zation directions of both electrons and holes are chosen 
along the z direction. It is obvious that the phonon states 
play a key role in these transitions. For transitions from 
HH and LH bands, the LA and LO phonon-assisted pro- 
cesses inject spin polarization along the —z and z di- 
rections, respectively, with a DSP of 50%, yet there is 
no spin polarization from the TA/TO phonon-assisted 
processes. Despite its spin independence, the electron- 
phonon interaction still affects the selection rules for spin 
injection since all states involved are not pure spin states 
due to spin-orbit coupling. Therefore, it is not adequate 
to treat the indirect transitions as a spin dependent vir- 
tual optical transition combined with a phonon emission 
or absorption process that does not affect the spinj 26 i 27 

Now we look at the transitions to the electron states in 
the X valley. The possible optical transitions are compli- 
cated, and are shown in Fig. [31 To simplify the diagram, 
we choose the quantization axis of the hole states along 
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FIG. 2: (color online). Diagram for optical indirect transi- 
tions from the band edge of the valence bands (HH, LO, and 
SO) to the band edge of the conduction band (CB) in the Z 
valley, under a~ light. Both the hole and electron states are 
quantized along the z direction. The red arrows stand for al- 
lowed transitions, whose probabilities are proportional to the 
product of the numbers next to the arrows and factors of |lf , 

§T 4 2 , and |[T 2 (/) ] 2 for LA, LO, and TA/TO phonon-assisted 
processes, respectively. 



the x direction, and that of the electron states along the 
z direction. Here the LO and TA/TO phonons can play 
a role in spin injection, but the LA phonon cannot. Thus 
there is strong valley anisotropy in the injection of spins. 
In Li and Dery's approximations, the TA/TO phonon- 
assisted processes give a DSP of 1/3 from the HH and 
LH bands, and the LO phonon-assisted process gives no 
spin polarization. In our more detailed analysis, we can- 
not simply identify a DSP for each process, as some of 
these processes are determined by more than one nonzero 
parameter. We discuss the actual values of these at the 
end of section IIVI 

Because the band edge transitions strongly depend on 
the choice of the electron and hole states, it is construc- 
tive to give the nonzero components of A a ^. CVT . They can 

be specified by giving the values of A C vt , in terms of 
which A a z. CVT are found following the pattern of Eq. J8|) 
and © . The results are shown in Table |IIJ In the cal- 
culation, we have used the result that the contributions 
for the HH and LH bands at the T point are exactly the 
same. 



From Table HU we find that some results are similar to 
those for direct gap injection: The band edge transition 
magnitudes for carrier injection are the same for each 
valence band, while the ones for spin injection are only 
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LO 


r i 

(±3/2)' \ 2 / (±1/2)' 


t 1 

IX 

(±3/2)' \ (±1/2)' 
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(±172)' 
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t 


1 


(±372)' \ / ( ±my 
\ ; 

so 
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FIG. 3: (color online). Diagram for optical indirect tran- 
sitions from the band edge of the valence bands (HH, LO, 
and SO) to the band edge of the conduction band (CB) 
in the X valley, under a~ light. The electron states are 
quantized along the z direction, and the prime beside the 
hole states indicate that they are quantized along the x di- 
rection. The transition strength for each color arrows in 
each phonon branch is as follows: In LA-assisted transitions, 
the factor is for black arrows. In LO-assisted transi- 

tions, the factors are §T 4 for black arrows, i(T 4 ±2T 5 ) 2 for 
blue thin dotted (+) and red thin dashed (— ) arrows, and 
Yj(T4 iTs) 2 for brown thick dashed (+) and green thin dot- 
ted (— ) arrows. In TA/TO-assisted transitions, the factors 
are \[t[' ] ] 2 for black arrows, i [{T^f + 2(T 2 (,) ) 2 ± 2T 2 (,) T 1 (,) ] 
for red thin dashed (+) and blue thin dotted (— ) arrows, and 
^[2(T^) 2 + (T 2 (,) ) 2 ± 2T 2 ( ' ) T 1 (,) ] for green thick dotted (+) 
and brown thick dashed (— ) arrows. 



T 


TA/TO 


LA 


LO 


c 


so 


(0,-2TWT 2 W) 


(-2T 3 2 ,0) 


(2T 4 2 , 2t 4 T 5 ) 


LH,HH 


(o,±tW) 


(|T|,0) 


/ 1 rri2 1 rji rj-i \ 

\ T\ 1 i ' ~ a i 4-t5) 


e 


SO 


(2[T 2 <"] 2 ,i[7f] 2 ) 


(|T|,0) 


(2T 4 2 , 2t 5 2 ) 


LH,HH 


( |[r 2 W]M[TW] 2 ) 


(|T 3 2 ,0) 


(|T 4 2 , §T 2 ) 



TABLE II: Band edge transition rates (A^} T , A^ivr) with A 
being £ for carrier injection and ( for spin injection. 



identical for HH and LH bands, and satisfy 

ESv 2) -0. (18) 

V 

The vanishing sum does not mean that the spin polariza- 
tion becomes zero when the laser pulse is wide enough to 
involve all valence bands, because the densities of states 
for different valence bands are different. However, as 
the photon energy increases, the optical transition oc- 
curs away from the band edge, where the involved elec- 
tron and hole states are superpositions of the band edge 
states of all valence bands. As a consequence of Eq. (|18l) . 
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important band mixing leads to a small DSP. 



IV. MODEL FOR ELECTRON STATES AND 
PHONON STATES 

To look at the carrier and spin injection away from the 
band edge, a full band structure model of the electron and 
phonon states is necessary. We use the EPM22r— for elec- 
tron states and the ABCM 32 for phonon states. We de- 
scribe these now, and give the resulting electron-phonon 
interaction for indirect gap injection. In the EPM, elec- 
trons are described by the pseudo-Hamiltonian, 



P 

2mo 



^v(r - R ia ) 



(19) 



Here p is the momentum operator; Ri a = Ri + t q is 
the equilibrium position for the a th (a = 1,2) atom 



the 



primitive cell located at Ri with t\ = 



and T2 = f (1,1,1), and a is the lattice constant; 
v(r,p) = v L (r) + vu L (r,p) + v so (r,p) is the atomic em- 
pirical pseudopotential, in which i'i(r) and v NL (r,p) are 
the local and nonlocal spin independent pseudopoten- 
tials, respectively ) 29 ' 30 and v so (r,p) is a spin dependent 
contribution 31 introduced to fit the spin split-off energy 
at the r point. 34 Eigenstates of H e are given by the 
Bloch states \nk) = J2 g c nk{g)\k + g)i where |fc) arc 
plane wave states, g is a reciprocal lattice vector, and 

c-nk — fc£ fe c^ lfc ^ are coefficients obtained by solving 
the single-particle Schrodinger equation. 

The Fourier transform of the pseudopotential v(r,p) 



v(k 1 ,k 2 ) 



dr 



—iki-r 



v(r,p)e 



ih 2 -r 



(20) 



For the local pseudopotential Vi(r) taken to be of the 
form Vi,(r) with r — \r\, we have that i>i,(fci) k 2 ) depends 
only on |fei — ka|, and we write it as v L (\ki — A^l). Chc- 
likowsky and Cohe n 29 ' 30 showed that a suitable electron 
band structure can be produced by taking the values of 
v L (k) only at k 2 = 3, 4, and 11 (2ir/a) 2 , which lead to 
the conduction band energy at the X point of 1.17 eV. 
In order to produce the direct band gap E g , the indirect 
band gap Ei g , and the location of the conduction band 
edge correctly after including the spin-orbit coupling, 
we use the following parameters Vi(k) = —3.496, —0.544, 
0.437, 0.429, and 0.1373 eV for k 2 = 3, 4, 11, 16, and 
19 (27r/et) 2 . The calculated band structure is shown in 
Fig. [TJ The calculation is for zero temperature. With 



increasing temperature, the electron-phonon interaction 
induces changes at the band edge and shifts the indirect 
and direct band gaps.— ~— For silicon, the shift in the 
conduction band edge is minor; we absorb it into Ei g in 
the following. 

For the phonons, the ABCM is used for calculating 
the polarization vectors e qX and energies Ml q \. The cal- 
culated band structure fits the experiments welli^ 2 - The 
energies of phonons with wave vector which are in- 
volved in the band edge transition, are hfl^o TA = 19 meV, 
f&lfco la = 43 meV, h£l k o LO = 53 meV , and hQ k o T0 = 
57 meV. 

By shifting the atom position from the equilibrium po- 
sition Ri a to Ri a + Ui a , and then expanding the elec- 
tron Hamiltonian (|19[) to linear order in Ui a , we identify 
the electron-phonon interaction as H ep = — J2ia ^ v { r ~ 
Ria) ■ Uia- The atomic displacement is usually expanded 
by the phonon polarization vectors as 



qel st BZ,A 



1/2 



{a qX + al qX )e% x e^ , (21) 



with p being the mass density of silicon. The transi- 
tion matrix elements between different electron states are 
given as 



9192 



A — iAk- 

q,ot 



xc l 2 k 2 (9i)v(ki fc 2 + g2)c n2 k 2 (g2) , (22) 

where q = k\ — k%, Ak = k\ + g\ — k-2 — 92- 

In calculating the matrix elements of the electron- 
phonon interaction given in Eq. (|22p . the values of v^(k) 
at all k points are necessary, and they are obtained us- 
ing the natural cubic spline interpolation 4 ^ on the dis- 
crete points given above and two further restrictions: 
the first is the value i>l(0); the second is a cut off at 
high k, v L (k > k cut ) — 0, where k cut is a cut-off value. 
Bednarek and Rossler 20 showed that the values assumed 
for «l(0) and k cut strongly affect the calculated matrix 
elements of the electron-phonon interaction in Eq. (|22l) . 
In our calculation, we find that the assumed value of 
k cut doesn't significantly affect H ep for k cut > 3kp, and 
we set kcut — 3fcp. We show the dependence of indi- 
rect absorption on v L (0) in Tabic Hill and HV1 by taking 
«l(0) = — §-Ep (denoted as case A) and ^l(0) = (case 
B). Here /cf and Ep are the Fermi wave vector and Fermi 
energy of the free electron gas appropriate to the valence 
electron density in silicon. 



For carrier injection, the contributions from the 



phonon branches other than TA are essentially the same 
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TA 


TO 


LA 


LO 




If 


TiT 2 


[T{f 


(U) 2 


Tin 


Ti 


Ti 


Ti 


T 4 T 5 


A 


0.020 


0.018 


-0.019 


0.131 


1.000 


-0.358 


0.018 


0.218 


0.001 


-0.014 


B 


0.066 


0.090 


-0.077 


0.120 


1.024 


-0.348 


0.020 


0.192 


0.026 


-0.070 



TABLE III: Relative band edge values of all quantities listed in Table ILT1 for the phonon emission process. The photon energies 
for different phonon branches are taken as frw = Ei g + hO, k a x . All values are normalized with respect to the value of |T 2 | 2 in 
case A. Its value is (Ti) 2 = 1.8 x 10~ 86 J 2 V~ 2 m _:L . 





Z valley 


X valley 


total 


c z = 

?cHHA + 


s z = 

;,, l >cHHA+ 


DSP Z = 

s z /c z 


c x = 

2 VScHHA4- ' ?cHHA + / 


s x = 

2 7(2) 

h ^cHHA+ 


DSP X EE 

8 X /C X 


C EE 

4c x + 2c z 


S EE 

As x + 2s z 


DSP EE 

s/c 


A 


TA 


0.0117 








0.019 


— 0.006 


-32% 


0.100 


-0.025 


-25% 


LA 


0.0122 


-0.0062 


-51% 


0.006 


0.00022 


4% 


0.0488 


-0.0115 


-24% 


LO 


0.145 


0.0727 


50% 


0.073 


0.00473 


6% 


0.583 


0.164 


28% 


TO 


0.667 


-0.0056 


-1% 


0.42 


-0.119 


-28% 


3.01 


-0.489 


-16% 


B 


TA 


0.0603 


-0.00045 


-1% 


0.074 


-0.026 


-35% 


0.42 


-0.103 


-25% 


LA 


0.0135 


-0.0067 


-50% 


0.0068 








0.054 


-0.013 


-24% 


LO 


0.128 


0.064 


50% 


0.072 


0.023 


32% 


0.546 


0.222 


41% 


TO 


0.682 








0.421 


-0.116 


-28% 


3.05 


-0.475 


-16% 



TABLE IV: Band edge values of ^hha- The classification in Z and X valleys follows Eq. (|23p . All values are normalized with 
respect to the value of IT^] 2 in case A. 



in case A and B, while for spin injection the phonon pro- 
cesses involving TA or LO phonons and electron injec- 
tion into the X valley give contributions that depends 
strongly on the value of i>l(0). The DSP for the LO 
phonon-assisted process is about 6% in case A and 32% 
in case B. Nevertheless, by far the most important con- 
tribution to indirect gap injection comes from the TO 
phonon-assisted process, which is almost independent of 
Wl(0) for both carrier and spin injection. In the following, 
we use the parameters of case A. 

Contrary to what was assumed in Li and Dery's re- 
sults, our results show that T2 = 2T\ or T' 2 — 2T[ cannot 
be well satisfied. But the calculated DSP for TA/TO 
phonon-assisted processes in the X valley is still close to 
1/3. The value of T5 is reasonably small. 



V. CALCULATIONS 

To obtain the injection rates, we can focus only on 
energies near the band edge and rely on Eq. (fTTj) . or 
numerically evaluate Eq. (jSJ) using the results of a full 
band structure calculation. Both strategies require a 
six-dimensional integration over the electron and hole 
wavevectors ranging within the BZ. In the present work, 
we vary the photon energy fouj to about 1.5 eV above Ei g , 
which results in an effective integration volume consist- 
ing of about 3/8 of the volume of the whole BZ, and a 
very demanding calculation. Similar to the integration 
used in the calculation of direct gap injection, the inte- 
grand here is composed of an energy conservation term, 



i.e. the Dirac 8 function, and a transition matrix ele- 
ment term. In a calculation of direct gap injection, the 
5 function is evaluated using a LAT M 24 i 35 on a very fine 
grid, and the transition matrix elements are calculated 
on each grid point. However, this method is not prac- 
tical in a calculation of indirect gap injection, because 
the transition matrix elements are too complicated to be 
evaluated on each point of a fine grid. Instead, we find 
that we can obtain a converged result by using separate 
grids for these two terms, adopting a finer grid for the 
Dirac S function and interpolating the transition matrix 
elements on a rougher grid. The details of this method 
are in Appendix |XJ 

In our calculation, the valence bands include HH, LH, 
and SO bands, the conduction bands include the lowest 
and the first excited conduction bands, and the interme- 
diate states are chosen as the lowest 30 bands to ensure 
convergence. 

The total carrier and spin injection from any pair of 
valence bands via an emission or absorption process in- 
volving a phonon of any branch, as well as the carrier and 
spin injection into any conduction band valley, are com- 
pletely determined once the quantities £^a± an( ^ ^cv\± 
are specified: The nonvanishing Cartesian tensor compo- 
nents follow from Eqs. ([5]) and ©, the sum over the dif- 
ferent valleys follows from Eq. (ITOl) , and the full response 
tensors follow from Eq. @; once these are determined 
the injection rates can be calculated from Eq. ([T]) for any 
light polarization. Instead of presenting our calculated 
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(1 2) (1 21 

results for Q v \± and Qv\± bekwj we present instead 



AX) 

$cv\± 



^cvX± 

Ax) 

^cvX± 

AZ) 

> >cv\± 



AD 

Sct)A± 



A2) 1 
' Vi,A±J 



AD 

^cvX± ' 

_A2) 
^cvX± 

^cvX± 



(23) 



These are the terms that appear in a simple excitation 
scenario using o~ light, as described in the following sec- 
tion. They correspond to the injection into different val- 
leys via the different processes. Nonetheless, we stress 
that given the quantities in Eq. (|23[) . we can construct 
the carrier and spin injection in Eq. ([1]) for any polariza- 
tion using Eqs. ©, ©, and ©. 



VI. RESULTS 

In the following, we focus on injection under o~ light. 
In this case, the six valleys can be divided into two sets: 
{Z,Z} and {X,X,Y,Y}. The injection is identical for 
all the valleys within each set. The carrier and spin in- 
jection from the valence band v into the conduction band 
c in the I th valley via a A^-branch phonon emission (+) 
or absorption (— ) are identified by fij- cv \± and Si- cv x±, 
respectively, which are given from Eq. (|2"3"|) as 



ni-cv\± 

°I;cvX± 



in 



(Z'xJEol 



(24) 



with Sj!c V x± = 0- Accordingly, the coefficients 



I;cvX± 
cu±;AEr 



ad = v 

cnA±' Scut / t ^cvX± 



±;AGr 



are also used in the following to understand the injection 
properties. Here r = TA, LA, LO, TO are the phonon 
branches; as before, X)agt indicates to a sum over all 
phonon modes in the r th branch. Similar notation is 
used for the injection coefficients of spins. We find that 
the polarization direction of injected spins in each valley 
is parallel or anti-parallel to the z direction. The DSP is 
defined as 



DSP = 



fi£/2' 



DSP 



(i) 
{-} 



AD 



(25) 



Here {•••} indicates the subscripts for DSP^) are the 
same as that for £W and 



A. Carrier injection 

1. Band edge carrier injection at 4 K 
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FIG. 4: (color online). Spectra of carrier injection rates £ 
(solid curves), £^ z ' (dashed curves), and £^ x ' (dot-dashed 
curves) at 4 K. The black thick and red thin curves are the 
results calculated from Eq. ([6]) and Eq. respectively. 

We first study the carrier injection at the band edge 
at 4 K. For band edge injection, only the transitions be- 
tween the lowest conduction band and the valence bands 
need to be considered; the first excited conduction band 
is ignored due to its small density of states. At this low 
temperature, only the phonon emission process is impor- 
tant. In Fig. |4] we show the spectra of the total injection 
rate (black thick solid curve) as well as the injection rates 
in Z and X valleys (black thick dashed and dot-dashed 
curves). All injection rates increase with increasing pho- 
ton energy. These results are consistent with the an- 
alytical results in section IIII1 where the injection rates 
around the band edge are approximately proportional to 
the JDOS. The difference between the injection rates in 
the Z and X valleys shows that the injection is valley 
anisotropic. To understand the contribution from each 
phonon branch, we plot phonon-resolved injection rates 
in Fig. [5j We find that each phonon-resolved spectrum 
has a shape similar to the total. In our calculation, the 
importance of the phonon-assisted processes are in order 
of TO > LO > TA > LA, in which the contribution of 
the LA-assisted process is less than 5%. 

At the band edge, the injection rates can also be ob- 
tained from the simplified calculation given by Eq. (jlll) , 
which can be used to identify how well the selection 
rules in Table |H] work. The results of the simplified cal- 
culation are also plotted in Figs. 0] and [5] as red thin 
curves. Compared to the full calculation, Fig. [5] shows 
that these results have smaller values for the TA, LA, 
and LO branches, and a larger value for the TO branch. 
The differences arc significant for acoustic phonons, and 
minor for optical phonons. Even for optical phonons, 
the difference is about 30% at excess photon energy of 
0.25 eV. However, because the errors for the two most 
important phonon branches, i.e. TO and LO branches, 



9 



25 0.3 




0.05 0.1 0.15 0.2 0.25 0.3 
hu - E, t (eV) 



0.05 0.1 0.15 0.2 0.25 0.3 
&j - E t (eV) 



FIG. 5: (color online). Phonon-resolved injection rates £f 
(solid curves) and ff (dashed curves). The black thick and 
red thin curves are the results calculated from Eq. ((6} and 
Eq. (fTT|) , respectively. 



mined only by the density of states of each valence band. 
Similar results are obtained for the other three phonon 
branches. 



2. Comparison with experiment at 4 K 




1.16 1.2 1.24 1.2c 
hu (eV) 



are opposite, the difference in the total injection rates be- 
tween these two methods is not so great. This difference 
illustrates not only the variation of the transition matrix 
elements on k c and k v , but also the failure of the sim- 
plified formula at high photon energy. Nevertheless, the 
simplified formula gives the correct qualitative results. 



FIG. 7: Spectra of y aftw. The solid curve is our theoretical 
result. The red thin curve is experimental data. 5 




0.05 0.1 0.15 0.2 0.25 0.3 
hw - E, g (eV) 

FIG. 6: (color online). The contribution from each valence 
band to the injection rates £' z ' (thick curves) and £ x ' (thin 
colored curves) for the TO assisted process. 

To better understand the details of carrier injection, 
we also plot the contribution from each valence band 
to ^ x ' z) for the TO-assisted process in Fig. [6j There 
are two important features: First, all injection rates in- 
crease with photon energy, which can be understood by 
the increase of the involved JDOS with photon energy. 
Second, the HH band gives the greatest contribution to 
the injection rate, and the SO band gives the smallest. 
As the transition matrix elements are the same for these 
three bands at the band edge, the magnitude is deter- 



Fig. Ogives the photon energy dependence of y/a(oj)huj 
at T — 4 K. Here a(ui) = Sw^(w)/(2n(w)ce ) is the ab- 
sorption coefficient, n(uj) is the refractive index of silicon, 
c is the speed of light, and eo is the vacuum permittiv- 
ity. Our result (the solid curve) shows two important 
features of the indirect gap injection: Near the onset of 
absorption, for hui < Ei g + hfl^ L0 , the TA phonon emis- 
sion process dominates, then the optical phonons make 
important contributions for fku > Ei g + fiwfco LD . The 
separation between these two regions is indicated by the 
kink in the figure. In each region, the lineshape of v 'ahui 
is approximately a linear function of Huj, which is consis- 
tent with the results in the parabolic band approxima- 
tion, i.e., \J ahaj oc fkj — Ei g — Klfj3\ i 11 - ' 15 With respect to 
these features, our results match the experimental results 
very well. 

Yet compared to the experiment by Macfarlane et a/.,— 
our result fails to show the correct energy shift and line- 
shape at the beginning of each region. Both of these fea- 
tures are related to the excitonic effect, which is absent in 
our calculation. When the excitonic effect is considered, 
both the bound and continuum exciton states contribute 
to the absorption, as discussed in detail by Elliott ^ The 
onset of absorption is shifted to lower energy by the pres- 
ence of exciton bound states, and modified to give a line 

shape a oc ^jhw — (Ei g + h£lf.o\ — E ex ) associated with 

the first exciton bound state and each phonon branch, 
with an exciton bind energy E ex sa — 14 meV^ the line- 
shape of the absorption associated with the exciton con- 
tinuum states is similar to that calculated without in- 
cluding the electron-hole interaction. Despite all this, in 
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the next we will see that the absorption at high photon 
energy can still be well described by our model. Even 
at low photon energy, we expect our model can give a 
reasonable description of the DSP, because the excitonic 
effects are to good approximation spin independent. 



3. Carrier injection at high temperature 
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FIG. 8: (color online). Spectra of £, and at 77 K 

(black thick curves) and 300 K (red thin curves). 

The indirect gap injection rates depend on the phonon 
number, and thus on the temperature. With increasing 
temperature, the phonon absorption processes come into 
play, and the onset of the injection spectrum moves to 
lower photon energy. The equilibrium phonon numbers 
N k o x at different temperatures are as follows: At 77 K, 
all phonon numbers can approximately be ignored, so 
only the emission processes occur. At 300 K, the phonon 
numbers become 0.92 for TA, 0.23 for LA, 0.15 for LO, 
and 0.12 for TO. Because the TA phonon has smaller en- 
ergy, the absorption induced by the TA phonon-assisted 
process is the most temperature-sensitive. 

In Fig.[8]the spectra of £, and are plotted for 
77 K (black thick curves) and 300 K (red thin curves). 
All injection rates increase remarkably with temperature: 
£ displays about 50% increase when the excess photon 
energy is 1.5 eV. Note that the injection rates at 77 K and 
4 K are approximately the same. The phonon-resolved 
absorption spectra at different temperatures are plotted 
in Fig. [SJ The injection rate from each phonon branch 
increases with temperature, and the increment is most 
significant for the TA phonon branch. At 300 K, is 
comparable to Ctc? > a ^ ^ ow temperature, is l ess than 
half of 4o } - 
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FIG. 9: Phonon-resolved £ (z) and at 77 K and 300 K. 
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FIG. 10: (color online). The comparison of a between our 
theory (thick curves) and experiments^, (thin colored curves) 
at 77 K (solid curves) and 300 K (dash-dotted curves). 



Compared to the results near the band edge, the calcu- 
lations without including excitonic effects fit the experi- 
ments better at high photon energy. The difference be- 
tween theory and experiment at low photon energy can 
be attributed to excitonic effects and the temperature 
dependence of the indirect band gap energy, neglected 
in our calculation which takes that energy as its zero 
temperature value. Clearly our neglect of excitonic ef- 
fects and of this gap energy shift has less consequence at 
higher energies. In the calculation, we use the experimen- 
tal frequency-dependent refractive index n(ui)§- Because 
of the increase in phonon number, the absorption coeffi- 
cient at 300 K is remarkably larger than at 77 K. 



5. Temperature dependence of carrier injection 



4- Comparison with experiment at high temperature 

In Fig. 1101 we compare our results for the absorption 
coefficient a with experimental results^ at 77 and 300 K. 



Fig. [TT] (a) gives the temperature dependence of the ra- 
tio£(T)/£(0) (black solid curves), ^ X) (T)/^ X \0) (green 
dash-dotted curves), and ^ (T)/£ (z) (0) (red dashed 
curves) for hu - E ig = 70, 100, and 1000 meV. The 
injection rates increase with temperature for all photon 
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FIG. 11: (color online), (a) Temperature dependence of 
£(T)/£(T = 0), e X) (T)/e X) (0), and (T)/£ (Z) (0) for ex- 
cess photon energy at 70, 100, and 1000 meV. (b) Temper- 
ature dependence of phonon resolved ratio for excess photon 
energy at 100 meV. Black solid curves: the total injected car- 
riers, green dash-dotted curves: values in the X valley, red 
dashed curves: values in the Z valley. 
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FIG. 12: (color online), (a) Spectra of spin injection rates 
£ (solid curves), Q z ^ (dashed curves), and Q x ' (dash-dotted 
curves) at 4 K. (b) The corresponding DSP spectra. The la- 
belled energies are hoj2,3,4 — Ei g = 19, 57, 90 meV, respectively. 



B. Spin injection 



energies, and their slopes decrease with photon energy. 
When the temperature is lower than 70 K, the injection 
rates are almost independent of temperature, while they 
become nearly a linear function of temperature at tem- 
perature higher than 200 K. Their slopes depend on the 
photon energy with larger slopes at lower photon ener- 
gies. The slopes for the injection in the Z and X valleys 
are found to be different at low photon energy, with the 
difference tending to disappear at high photon energy. 
The temperature dependence of phonon-resolved injec- 
tion rates for excess photon energy 100 meV are plot- 
ted in Fig. QT] (b) . The temperature dependence of the 
phonon-resolved injection rates are similar to the total, 
but show different slopes for different branches. In order 
to give an comprehensive understanding of these results, 
we use the simplified formula for the band edge injection 
in Eq. (jlip . At high temperature, the phonon number 
Nj^ox is approximately fe#T '/Mlk°X) an d Eq. (jlll) for ab- 
sorption near the band edge can be rewritten as 



= 1 



k B T 



Huj- E, 



hw - E ig - HQ k o X 



(26) 

The linear dependence at high temperature is obvious. 
The slope is determined by the photon energy, and de- 
creases as the photon energy increases. At photon en- 
ergy displayed in Fig. [TT] (b), however, the dependence 
of the slope on phonon energy is too complicated to be 
described by this simple formula. 



1. Band edge spin injection at 4 K 

In Fig. [12] (a) we show the spectra of the spin injec- 
tion rate in the Z valley (dashed curve) and the X valley 
(dash-dotted curve) as well as the total (solid curve) at 
4 K. The spin injection rates increase with photon en- 
ergy, and show strong valley anisotropy, as do in car- 
rier injection. But the photon energy dependence of the 
spin injection rates does not follow the quadratic JDOS, 
but nearly a linear function. This indicates the strong 
wavevector dependence of Cct c vk v \ gi yen i n Eq. flU). The 
injected spins in the X and Z valleys have opposite po- 
larization direction. Fine structures of spin injection can 
be found from the DSP spectra given in Fig.[T2](b). The 
total DSP at the band edge (fiw 2 — E ig = 19 meV, cor- 
responding to the TA phonon emission process) is about 
—25%, while the DSP in the Z and X valleys approxi- 
mates and —32%. When the photon energy is higher 
than hu>2, the spectra can be divided into three regions: 
i) Huj € [hu>2, Hui^\ with Hui^ — Ei g = 57 meV (= flf.o xo ). 
Here the DSP and the DSp( x ) decrease rapidly with in- 
creasing energy to a minimum value of —1% and —5%, 
respectively, while the DSP^' increases to a maximum 
value of 5% when the photon energy is slightly smaller 
than hu>3 , and then decreases slightly to a local minimum 
at Hu>3. ii) Hw € [thu)3 , huj^] with Huj^ — Ei g = 90 meV. 
Here the DSP and the DSP (X) increase to a maximum 
value -4% and -10%, and the DSp( z > first slightly in- 
creases and then decreases, iii) Huj > H104. All DSP de- 
crease monotonically to zero. The special photon energy 
Huj 3 is related to the TO phonon energy given in Sec. IIV[ 
the features here are formed by the contributions from 
different phonon branches, which are plotted in Fig. [T51 

According to the phonon energies, the injection edge 
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FIG. 13: (color online). Phonon-resolved spin injection rates 
[(a) and (c)] and DSP [(b) and (d)] in the Z [(a) and (b)] and 
X [(c) and (d)] valleys. 



FIG. 14: (color online). Phonon and valence band resolved 
DSP in the Z valley (thick curves) and X valley (thin colored 
curves) . 



energies of different phonon branches are ordered as TA 
< LA < LO <TO, which is obviously shown in Figs. [T3l 
(b) and (d). For the spin injection in the Z valley given 
in Fig. [13] (a), the contributions from all phonon branches 
are almost of the same order of magnitude, even including 
the LA phonon branch, which gives a negligible contribu- 
tion to carrier injection. Among these branches, the LO 
branch gives the largest contribution. The correspond- 
ing DSP are shown in Fig. [13] (b). The DSP induced 
by TA branches starts from zero, reaches a maximum at 
about 45 meV above the onset of absorption, and then 
decreases, while the LA phonon branch contributes a neg- 
ative spin injection rate. LO and TO phonons take effect 
for hjj > Hlj 3 . Therefore, for the total DSP in the Z 
valley (in Fig. I12[) . the first increase to the maximum 
is induced by TA phonons, the following dip is induced 
by LA phonons, and the second peak is induced by LO 
phonons. Figs. [131(c) and (d) give the spin injection rates 
and DSP in the X valley. Here the contribution from 
the TO phonon branch dominates the spin injection, and 
those from other phonon branches only give minor con- 
tributions. As in the Z valley, the spin injection in the 
X valley (in Fig. [T2|) starts from the TA phonon branch; 
its DSP decreases with photon energy from a nonzero 
band edge value, and gives the fast decrease for photon 
energy in hw G [tkJ2, TujJs]- At hjj 3 the TO phonons come 
into play, leading to injected electrons with a large spin 
polarization. As the contribution to the total spin po- 
larization from these electrons begins to dominate, the 
DSP reaches a maximum at Hui^, and decreases at high 
frequencies. 

We now turn to resolving the spin injection rates into 
contributions from different valence bands. The results 
are similar to the corresponding resolution of carrier in- 
jection rates. The HH band gives the largest contribution 
due to its large density of states, while the SO band gives 
the smallest. However, there are subtleties in the DSP 
of injected spins from each valence band. In Fig. 1141 we 



give the valence band and phonon-resolved spectra of the 
DSP. While the band edge values of the DSP are consis- 
tent with results in Table HVl the spin polarization direc- 
tion shows a complicated energy dependence based on 
whether or not the band edge DSP is zero. For processes 
with nonzero band edge values, the DSP decrease from 
the nonzero values. For processes with zero band edge 
values, the DSP increase to a maximum value, and then 
decrease. Some of these DSP change sign with increasing 
photon energy. 



2. Spin injection at high temperature 
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FIG. 15: (color online). Spectra of spin injection rates and 
DSP at 77 K and 300 K. hu)i - E ig = -57 meV. 

In Fig. [I5]we plot the spectra of spin injection rates at 
77 K and 300 K. In contrast to the carrier injection rates, 
which increase considerably from 77 K to 300 K, both 
the total spin injection rates and the rates in each valley 
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change little from 77 K [Fig. [15] (a)] to 300 K [Fig. [15] (c)]. 
In order to show clearly the spin injection at the injection 
edge, the DSP are plotted in (b) and (d) respectively, the 
phonon absorption process becomes increasingly impor- 
tant as the temperature rises, and is shown by the left 
shifts of the onset of DSP to Hu>i, which is determined by 
the TO/LO phonon absorption processes. In addition, 
the peak appearing at H0J4 and 4 K in Fig. [T2] (b) be- 
comes fairly obscure at high temperature. These results 
can be better understood from the phonon-resolved DSP 
spectra at 300 K, given in Fig. [TBI 
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FIG. 16: (color online). Phonon-resolved spectra of spin in- 
jection rates and DSP at 300 K. 
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FIG. 17: (color online). Spectra of DSP for TO phonon ab- 
sorption and emission processes in the Z and X valleys at 
300 K. 



is induced by the phonon absorption processes, and de- 
creases with increasing photon energy. When the pho- 
ton energy is higher than hu>3, the TO phonon emis- 
sion process comes into play. However, the spin injection 
rate for phonon emission is approximately proportional 
to 1 + N^o TO) while that for phonon absorption is pro- 
portional to N k TOi thus the former gives much larger 
injection rates than the latter. This results in a peak 
in the total DSP for the TO phonon-assisted process, at 
energies where phonon emission becomes important. 



For the spin injection rates in the Z valley, shown in 
Fig. [TB] (a), the contribution from each phonon branch is 
still of the same order of magnitude, but the TO phonon- 
assisted process dominates the injection at high photon 
energy. For the spin injection rates in the X valley, shown 
in Fig. [TBJ (c), the contribution from the TO phonon 
branch dominates for all photon energies. However, the 
DSP, plotted in Figs. [TB] (b) and (d), show finer struc- 
ture in the low photon energy region than the DSP at 
4 K. The fine structure is induced by the combined ef- 
fect of the phonon emission and phonon absorption pro- 
cesses, which take effect at 300 K with considerable equi- 
librium phonon numbers for all phonon branches. Ac- 
cording to the simplified Eq. (fTTj) . the injection rates 
for phonon absorption/emission processes differ from the 
JDOS, J cv \(hoj =F hflko\), at low photon energy. There- 
fore, spectra of DSP for these two processes have similar 
lineshapes, except for a phonon energy shift left or right. 
This conclusion is confirmed by our numerical results, 
shown in Fig. [T7] Here the DSP from TO phonon ab- 
sorption and emission processes are plotted for spins in 
both Z and X valleys at 300 K, and all the absorption 
curves are right shifted by 20,^°. to- The overlap is obvi- 
ous. Based on these results, we return to Fig. [TH 

We consider the DSP from the TO phonon-assisted 
process in the X valley, shown as solid curve in Fig. [TB] 
(d), is taken as an example. The spectrum in [hu>i, Huis] 



3. Temperature dependence of spin injection 
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FIG. 18: Temperature dependence of DSP* 2 - 1 (dot-dashed 
curves), DSP^- 1 (dashed curves), as well as the total DSP 
(solid curves) at photon energies fkJi-4. 

Fig. HU shows the temperature dependence of the DSP 
at the different photon energies hiui^^. For photon en- 
ergy hwi, all DSPs are temperature-independent: The 
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TO phonon absorption process is the only process at this 
photon energy, and the injection rates of carriers and 
spins are approximately proportional to Nya xo'-> the two 
rates therefore scale the same with temperature, leading 
to a temperature-independent DSP. We check the tem- 
perature dependence of the DSP for all phonon branches 
at photon energy ftw — Ei g = 100 eV in Fig. [19] All 
these DSP are independent of temperature. This con- 
clusion is valid for all photon energies considered to this 
point in this paper. However, the total DSP, shown in 
Figs. HH] (b), (c), and (d), are temperature dependent. 
At the indicated photon energies, more than one phonon 
branch makes a contribution. The injection rates for dif- 
ferent phonon branches have different temperature de- 
pendences due to the different phonon energies. The ra- 
tio between the total injection rates of carriers and spins, 
i.e. the DSP, thus becomes temperature-dependent. At 
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FIG. 19: Temperature dependence of the DSP for TO and 
LO phonon emission assisted processes. 

high temperature, N q \ oc T is approximately valid for all 
phonon branches, so the injection coefficients of £ and £ 
are approximately proportional to T and give a constant 
DSP as saturation value. This is shown in Fig. [TSJ 



VII. CONCLUSION 

In conclusion, we have performed a full band structure 
calculation to investigate the indirect optical carrier and 
spin injection in bulk silicon. The injection spectra for 
carriers, spins, and DSP in each valley are studied in de- 
tail at different temperatures. The indirect gap injection 
is dominated by the transition from the HH band to the 
lowest conduction band due to the large JDOS. When in- 
cident light is propagating along one principal axis, the 
injection shows strong valley anisotropy. The injection 
rates induced by each phonon-assisted process increase 
with temperature. For carrier injection, we find that in 
the Z valley the TO phonon-assisted process dominates 
up to 300 K; in the X valley, it only dominates at low 



temperature, while the injection rates induced by the TA 
phonon-assisted process increases to a comparable value 
at 300 K. The higher the photon energy, the weaker the 
temperature dependence. For spin injection, we find that 
injected spins in the Z and X valleys have opposite po- 
larization directions. In the Z valley, the LO phonon- 
assisted process dominates around the band edge, while 
the TO phonon-assisted process dominates at high pho- 
ton energy; in the X valley, the TO phonon-assisted pro- 
cess dominates for all photon energies. 

The calculated absorption coefficients are in good 
agreement with experiments at high photon energy. Ex- 
perience from direct gap absorption might lead one to 
believe this would not hold, since first principle studies 
of direct optical absorption in Si and GaAs^r— demon- 
strated that the electron-hole interaction plays an im- 
portant role even for high photon energy. Full band 
structure calculations of indirect gap absorption, includ- 
ing electron-hole interaction, have yet to be done, and 
the effect of these interaction on the indirect absorp- 
tion coefficient at high photon energies have yet to be 
established. In any event, since to first approximation 
the nearly spin-independent electron-hole interaction will 
modify spin and carrier injections in the same way, it is 
reasonable to expect that the DSP of injected electrons 
in exciton continuum states will be insensitive to this in- 
teraction, which we neglect. Future work to confirm this 
is clearly in order, but in the interim we feel our study 
constitutes a good first investigation. 

The DSP spectra excibit a rich variety of behaviours. 
At 4 K, the maximum DSP is about —25% for photon en- 
ergy Eig+ClifO^A; absorption at this temperature and en- 
ergy is dominant by TA phonon emission process. With 
increasing temperature, the phonon absorption processes 
become important, and the DSP for this photon energy 
decreases quickly. At 300 K, the maximum DSP appears 
at the photon energy Ei g — f2 fc o T q as a value —15%, 
which only comes from the TO phonon absorption as- 
sisted process, and is temperature-independent. The 
DSP in the X valley can reach a maximum value —32% 
at 4 K and —26% at 300 K; both are larger than the total 
value and the value in the Z valley. Compared to bulk 
silicon, it should therefore be more efficient to inject spin 
in a confined silicon structure, where conduction band 
valleys are splitted. 
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Appendix A: Modified Tetrahedron Summation 



Both Eq. (|6] 
following form 



and the JDOS in Eq. (12]) arc of the 



/ dk c I dk v 8(e kc - e kv + hil kc - kv - Hoj) 
Jy JV 



xX(k c ,k v ) 



(Al) 



Here the integration spaces V' c and VI are the irreducible 
BZ wedges, which are confined by the T, X, L, U(K), 
and W points of a fee Brillouin zone. Our numerical 
scheme is based on the improved LATM,— ^ We divide 
V' c {VI) into N c (N v ) small tetrahedra as V' c 



(VI = Ef=iK'j) with verti ces k cJl I = 1,---,M C 
(k v . j,J = 1, • • • ,M V ). By a simple transformation, 
Eq. (IA1[) becomes 



Y N ° V' ■ 

L^i—\ c,i 



J deXij(uj,e) , 



In the same way, we get w^^e) for Eq. (| A2|) . Then we 
find 



X l3 (uj,e) 



and 



E <iv 

nm=l 

xXfecii k v 



e) 



with 



AH = /de^W c>J (e)W^ J (?uj-M2 feM _ fe „ iiJ -e) 



jj 



xA(fc Cj /,fc„ i j) 



(A5) 



with 



Ay(a;, e) = / dk c S(e kc - e)Yj(k e , Tiw - h£l vj - e) 



and 



y,(fc c ,e) 



dk v 5(e kv - e)X(k c ,k v 



(A2) 



(A3) 



In obtaining the equations above, we approximate the 
phonon energy to be a fixed value hflij between the two 
small tetrahedra V' c i and V^j. Now the Eqs. (|A2|I and 
(|A3|) have the same shape, so we only discuss Eq. (|A3|) 
in the following. Though the tetrahedron method can be 
applied to Eq. (|A3I) directly, we avoid doing this because 
the transition matrix elements are obtained within the 
scheme of the EPM, which are very time-consuming in 
calculation for each (fc c ,/, feu.j) pair. However, the direct 
tetrahedron method needs more k points for convergence 
than is computationally feasible. Instead, since the de- 
pendence of the transition matrix elements on k c and k v 
is weak, we linearly interpolate A(fe c , k v ) in tetrahedron 



X(k C7 k v ) 



m— 1 



X(k c ,k™)F™(k v ) 



with k™j being the position of four vertices of this tetra- 
hedron, and F™j(k v ) being the linear interpolation func- 
tion related to the m th vertex. Then Eq. (|A3I) becomes 

4 

Y^e) = 5>™(e)A(fc c ,fc™.) , 



dk v 6{e vkv - e)F™(k v ) , (A4) 



W v ,j{e) 



E 



(A6) 



In the above equations, we replace ?if2fe« 4 _fem. by Mlij. 

In indirect absorption, the excited holes are located 
around the top of the valence band, while the electrons 
are around the band edge of the conduction band. So we 
use a different division for k c and k v . With the division 
points k c j and k Vi j, we use the Delaunay triangulation 
method from the CGAL package^ to set up the tetrahe- 
dra. In these tetrahedra, the weights w™j and 10™^ can be 
calculated. However, the integration includes a S func- 
tion, which is a fast-varying function in k space, and the 
present tetrahedron is too rough to obtain the weights 
with required precision. We refine this tetrahedron into 
smaller ones. 

In our calculation, the division depends on the band 
structure. For conduction bands and the heavy hole 
band, k c j or k Vy j are generated with their distance not 
larger than 0.05 x — , while for light hole and spin split- 
off bands, the distances become 0.04 x — and 0.02 x — , 

1 a a ' 

respectively. For the heavy hole grid, the accurate 
weights can be obtained by refining each tetrahedron as 
8 3 smaller ones (each edge is refined into 8 parts). 

The convergence of the results is examined for the 
injection from heavy hole band to conduction band by 
changing the division points distance to 0.025 x — . The 
difference between the injection rates from these two di- 
visions is less than 5%. 
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